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Abstract 

We extract the deep inelastic scattering cross-sections ratio R = ctl/o~t in the 
range 10 -4 < x < 10 _1 from F2 HERA data using very simple relations based 
on perturbative QCD. The result depends on only one parameter 5, being 
x~ s the behavior of the parton densities at low x, which has been determined 
recently with a good accuracy by the HI group. 



In recent years the behaviour of deep inelastic lepton-hadron scattering (DIS) at the 
small values of Bjorken variable x, has been intensively studied. The present letter is 
devoted to the behaviour of the ratio of cross-sections of the absorption of a longitudinal- 
and transverse-polarized photon by hadron: R = o~l/°~t, & t small values of x. The ratio R, 
which may be represented as the combination of the longitudinal Fl(x, Q 2 ) and transverse 
F 2 (x,Q 2 ) DIS structure functions (SF): 

is a very sensitive QCD characteristic because it is equal to zero in the parton model 
with spin— 1/2 partons (it is very large with spin— partons). However, modern DIS 
experimental data (see the review in ||) are not accurate enough to determine R(x,Q 2 ). 
In addition, at small values of x, R data are not yet available, as they require a rather 
cumbersome procedure (see ||, for example) for the extraction from the experiment. 

We study the behaviour of R(x, Q 2 ) at small values of x, using the HI data |3|, [|] 
and the method || of replacement of the Mellin convolution by ordinary products. By 
analogy with the case of the gluon distribution function (see [|, [7|, g]) it is possible to 
obtain the relation between F L (x, Q 2 ), F 2 (x, Q 2 ) and dF 2 (x, Q 2 )/dlnQ 2 at small x. Thus, 
the small x behaviour of the ratio R(x, Q 2 ) can be extracted directly from the measured 
values of F 2 (x, Q 2 ) and its derivative. These extracted values of R may be well considered 
as new small x "experimental data" c . Moreover, when experimental data for R at small 
x become available with a good accuracy, a violation of this exactly perturbative relation 
will be an indication of the importance of other effects as higher twist contribution and/or 
of non-perturbative QCD dynamics at small x. 

We follow the notation of our previous work in refs. f| [|. The singlet quark s(x, Qq) 
and gluon g(x, Qq) parton distribution functions (PDF) d at some Qq are parameterized 



by (see, for example, |T0 



p(x, Qq) = A p x (1 - x) Up (l + e p ^ + l P x) (hereafter p = s, g). (2) 

The value of 5 is a matter of controversy. The "conventional" choice is 5 = 0, which 
leads to a non-singular behaviour of the PDF when x — > 0. Another value, S ~ ~, was 
obtained in the studies performed in ref. JTTJ] as the sum of the leading powers of ln(l/x) 



in all orders of perturbation theory. Experimentally, recent NMC data [12] favor small 
values of 5. This result is also in agreement with present data for pp and pp total cross- 
sections (see ||13||) and corresponds to the model of Landshoff and Nachtmann pomeron 



14]j with the exchange of a pair of non-perturbative gluons, yielding 5 = 0.086. However, 
the new HERA data [|3], [| prefer 5 > 0.2. For example, the 5 value obtained recently 
by HI group seems to depend slowly on Q 2 values. Its average value increases from 
5 = 0.228 at Q 2 = 8.5 GeV 2 to 5 = 0.503 at Q 2 = 800 GeV 2 . 

Further, we restrict the analysis to the case of large 5 values (i.e. x~ s 3> 1 ) follow- 
ing recent HI data ||. The more complete analysis concerning to the extraction of the 
longitudinal SF Fl(x, Q 2 ), may be found in 0, where we took into account also the case 
5 ~ corresponding to the standard pomeron. 

Assuming the Regge-like behaviour for the gluon distribution and F 2 (x, Q 2 ) at x~ s ^> 1: 

g(x, Q 2 ) = x- 5 g{x, Q 2 ), F 2 (x, Q 2 ) = X - S s{x, Q 2 ), 
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we obtain the following equation for the Q 2 derivative of the SF F 2 e : 

dF 2 (x,Q 2 ) = 1 _ s ^ ^ (q) mQ 2 )+r S p{a) xmQ2) + 0{x 2 ) 



dlnQ 2 



F L (x, Q 2 ) = x~ s £ (rff(a) P(0, Q 2 ) + r 5 Lp (a) xp'(0, Q 2 ) + 0(x 2 



(3) 



p=s,g 



where r^ p (a) and r v Lp (ot) are the combinations of the anomalous dimensions of Wilson 
operators 7^ = a^ ri +a 2 ^ ri +0{a 3 ) and Wilson coefficients f aB% v (l+aR p £ v ) +0(a 3 ) 
and olB™ + 0(a 2 ) of the n "moment" (i.e., the corresponding variables extended from 
integer values of argument to non-integer ones): 



r L(") = "5 
e 



S.I) 

L 



l + a{R s £ v -B 3 ^ 



f 



aB™ 



+ 0(a 3 ). 
l + a(R 9 £ v -B^B^/B™ 



0{a 3 



rU<*) = <*y®» + « 2 (li^ + + 2/W 7 ) + 0(a 3 



(4) 



and 



./ 



+ 0(a 3 ), 



p'(0, Q 2 ) = —p(x, Q 2 ) atx = 0, 
ax 



where e = J2i e ] is the sum of squares of quark charges. 
With accuracy of 0(x 2 ~ s ), we have for Eq.(P|) 



dF 2 (x,Q 2 
dlnQ 2 



rl} S (Ur d 9(x/U, Q 2 ) + rlfF 2 (x, Q 2 ) + (r d ss - r^)x^~s\x, Q 2 ) 



,1+5 



1+S\l-S~ri 



+ 0(x 2 ' s ) 



(5) 



Fl(x, Q 2 ) = r^ g \i Lg )- 5 g{x/^ 9 , Q 2 ) + rl + s s F 2 (x, Q 2 ) + (r s Ls - r^x^s'ix, Q 2 ) 
+ 0(x 2 - 5 ), 



with £ 



,l+<5 /„5 



/r% and £ Lg = r\f /r{ g . 



sg ' sg 



From Eq. (H) and (|D one can obtain Fl as a function of F 2 and the derivative 



(6) 



F L (x,Q 2 ) = -e 



r 



i+<5 



dlnQ 2 



r 



1+5 



2-5 „ „.l-5\ 



(7) 



where the result is restricted to 0(x 2 s , ax 1 s ). To arrive to the above equation we have 
performed the substitution 



2 



and neglected the term ~ s'(x£, sg , Q 2 



This replacement is very useful. The NLO anomalous dimensions "fQ)' n are singular in 



sp 



both points, n = 1 and n = 0, and their presence into the arguments of p(x, Q 2 ) makes the 
numerical agreement between this approximate formula and the exact calculation worse 



(we have checked this point using some MRS sets |10[ of parton distributions). 

0.5 and 5 



Using NLO approximation of r]+ 5 and r]^ s for concrete values of 5 



0.3 



we obtain (for f=4 and MS scheme): 



if 5 = 0.5 

F L (x,Q 2 ) = 



0.87 



1 + 22.9a 



dF 2 (0.70x,Q 2 
dlnQ 2 



+ 4.17aF 2 (0.70x,Q 2 



+ 0(a 2 ,x 2 ' 



-\ax l ~ 



if 5 = 0.3 

Fl(x,Q 2 ) = 



0.84 



59.3a 



dF 2 (0A8x,Q 2 
dlnQ 2 



+ 3.59aF 2 (0.48a;,Q 2 ) 



0(a 2 x 2 - 



"fax 1 " 



(9) 



With the help of Eqs. ([!]) and fl8|)-(P) we have extracted the ratio R(x, Q 2 ) from HI 
1994 data |J, determining the slopes dF 2 /dlnQ 2 from straight line fits as in ref. [|J ||. In 
the present calculation only statistical errors have been taken into account, and we have 
used A|jg — 225MeV^ in the calculation of the running coupling constant a s (Q 2 ) at two 
loops. 

Figure la shows the extracted ratio R at Q 2 = 20 GeV 2 for two different values of 
the parameter 5. It also shows BCDMS [15] and preliminary CCFR (see fl~6| ) data points 



where the errors are very much larger. For comparison we have also plotted various 
predictions for R using QCD formulas at 0(a 2 ) [I7J and parton densities extracted from 



fits to HERA data. The large difference between the result from MRS(G) and the latest 
set MRS(Rl) |JT5|] shows, as it is expected, the large effect on R of the unknown of the 
gluon distribution at small x. 

In figure la one can also see that the result from MRS(Rl) fits very well the points 
obtained with 5 = 0.5 for the lowest x data, although it fails to account for the highest x 
bins. The calculation with MRS(D-) is also statistically compatible with our data. 

By other part recent theoretical predictions on R based on conventional NLO DGLAP 
evolution analysis of HERA data (LBY) |19| and on the dipole picture of BFKL dynamics 
(NPRW) pm , both finding values 5 = 0.3, lie closer to the data points obtained with 
5 = 0.3 Eq. (|). 

Finally Fig. lb shows R for 5 = 0.3 (the value favoured by HI data [RJ) and at three 
different Q 2 values in comparison with the SLAC R(1990) parametrization pl |. One can 



see the very good agreement at x < 10 -2 even if only the statistical errors are taken into 
account. 

Notice that the points at the same x and different Q 2 are correlated by the form in 
which the derivative term dF 2 / dlnQ 2 is determined. 

In summary, we have presented Eqs. ([!]) and (|7|)-(|9]) for the extraction of the ratio 
R = <jl/<jt at small x from the SF F 2 and its Q 2 derivative. These equations provide 
the possibility of the non-direct determination of R. This is important since the direct 
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extraction of R from experimental data is a cumbersome procedure (see 0). Moreover, the 
fulfillment of Eqs. (P, (0)-© by DIS experimental data is a cross-check of perturbative 
QCD at small values of x. Our formulas can also be used as a parametrization of R as a 
function of the most widely used phenomenological F 2 . 

We have found that the results depend on the concrete value of the slope S. In the 
case 5 = 0.3, which is very close to the values obtained by HI group |3[ at the considered 
Q 2 interval, we found very good agreement with the SLAC parametrization [^ij and also 
a relatively good agreement with the studies based on NLO DGLAP and BFKL dynamics 
(see |19| and p0| , respectively). However the calculation performed with the latest sets 



of HERA parton densities using perturbative QCD at second order (see MRS(Rl) curve 
in Fig. la) predicts an slightly higher value for R. 
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Figure captions 



Figure 1: The ratio R = olI°~t at small x. The points were extracted from Eqs. ([!]), 
(|) and (g) using HI |, [| data. The dashed-dotted line (NPRW) is the prediction of 
Saclay group |Mj based on the dipole picture of BFKL dynamics. The band represent the 
uncertainty from the DGLAP analysis of HERA data by [19]. It is also shown BCDMS 
data |15| points at high x and the preliminary CCFR data point from (T^j . The solid lines 
in Fig. lb are the SLAC R1990 parametrization at Q 2 = 8.5, 20 and 35 GeV 2 (lower 
curve corresponds to lower Q 2 value). 
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R(x,Q 2 ) 



R(x,Q 2 =20 GeV 2 ) 




